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We study the properties of a quantum critical point which develops in a BCS superconductor when 
pair-breaking suppresses the transition temperature to zero. The pair fluctuations are characterized 
by a dynamical critical exponent z = 2. Except for very low temperatures, anomalous contribution 
to the conductivity is proportional to vT in three dimensions, but to 1/T in two dimensions. At 
lowest temperatures, the conductivity correction varies as T 1//4 in three dimensions, and as In (1/T) 
in two. 
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PACS numbers: 71.10.Hf, 71.27.+a, 74, 74.20.Fg 



The possibility of quantum critical behavior in itiner- 
ant magnets has attracted great attention in recent years, 
in part, because quantum criticality affords the possibiL 
ity of a controlled study of non-Fermi liquid behavior.EJ 
At a quantum critical point, order parameter fluctua- 
tions develop an infinite correlation range in both space 
and timcEIB The coupling between these fluctuations and 
conducting electron fluid is able, under certain circum- 
stances, to eliminate the formation of well-defined quasi- 
particles in the electron liquid, giving rise to a new kind 
of metallic behavior. 

In this paper we discuss the possibility of quantum 
critical behavior in superconductors. A quantum criti- 
cal point implies a finite value of the electron interac- 
tion strengths. At first sight, this would appear to rule 
out the possibility of a superconducting quantum crit- 
ical point, for conventional superconductivity develops 
for an arbitrarily small pair interaction. If the transition 
temperature is driven to zero in a pure BCS supercon- 
ductor, the pairing interaction and the pair fluctuations 
are completely eliminated. Fortunately, this is not the 
case in the presence of pair-breaking, which cuts off the 
logarithmic singularity in the pair susceptibility, requir- 
ing that the pairing interaction reach a critical strength 
before superconductivity develops. In this paper we char- 
acterize the quantum critical behavior which develops at 
this special point. Our results can be tested experimen- 
tally on conventional superconductors, such as Ce-doped 
Las They may also provide a useful diagnostic tool for 
the understanding of unconventional, e.g. heavy fermion 
superconductors. 

We begin with writing down the effective action of the 
pair field A in the vicinity of a quantum critical point: 
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FIG. 1. Schematic phase diagram in three (D = 3) and in 
two (D — 2) spatial dimensions. Shaded regions correspond 
to the superconducting (SC) phase. The normal phase corre- 
sponds to the Fermi liquid (FL) in three dimensions, and to 
a two-dimensional disordered metal (2DM) in two. 



which is valid for T,vq, \v\ << To- Here T is tempera- 
ture, N(0) is the density of states at the Fermi surface, 
v is a Matsubara frequency, q is momentum, vf is the 
Fermi velocity, and Sx is the deviation of the control pa- 
rameter x from its critical value x c , where the transition 
temperature turns into zero. 

We obtained this, action by repeating the Abrikosov- 
Gorkov calculations for an s-wave BCS superconductor 
doped witlumagnetic impurities. Apart from the original 
calculations, the frequency and momentum dependence 
of the disorder-averaged pair propagator was kept and 
at the end both the pair field and the momentum were 
rescaled to give (l)-(2). 

The quantity To is the only characteristic energy scale 
of the effective action (l)-(2). In a BCS superconductor 
it is of the order of the pair-breaking rate, which at the 
critical point is of the order of the transition temperature 
in the clean system. Both quantities are much less than 
the Fermi energy e^. However, in a strongly correlated 
material, To may, in principle, be of the order of tp. 

These two limiting cases correspond to different 
physics. In the BCS case, the quartic term in (1) can 
be neglected at all experimentally relevant temperatures 
T well below To. By contrast, if T) ~ ep, the feedback 
of the quartic term dominates at T << Tq. In the latter 
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TABLE I. Order of magnitude and leading temperature 
dependences of the imaginary part of the electron self-energy 
£(u;;p = pf;T = 0) due to scattering by the fluctuations of 
A, the specific heat coefficient correction SC(T)/T and the 
conductivity correction Sa(T) in the weak coupling (BCS) 
limit (To << £f) and in the strong coupling limit (To ~ £f). 
The value of 00 corresponds to the residual normal state con- 
ductivity. 
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case, one has to regard To a s a phenomenological parame- 
ter, resulting from a strong coupling, or non-BCS pairing 
mechanism. 

With (l)-(2) at hand, one can calculate various ther- 
modynamic and transport properties of interest. Table I 
presents the results for the zero-temperature quasiparti- 
cle decay rate due to scattering by the A field, and for the 
leading corrections to low-temperature thermodynamics 
and transport at T c = 0. In the BCS limit, corrections 
to the specific heat coefficient come from Gaussian fluc- 
tuations of the pair field.tl The correction for a strong 
coupling case was found ind using renormalization group 
methods. The quasiparticle decay rate due to scatter- 
ing off pair fluctuations is estimated by the diagram on 
Fig. 2 (a). It is essential for the calculation of conductivity 
and of the quasiparticle decay rate, that the electron ver- 
tex corrections are not singular, since the pair-breaking 
makes the lifetime of a Cooper pair finite. The lead- 
ing conductivity., correction is given by the Aslamazov- 
Larkin diagramQ shown on Fig. 2(b), which can be re- 
garded as conductivity of particles with the inverse prop- 
agator given by (2), at Sx = 0. 

The calculation was done by renormalization group 
analysis of the expression for the Aslamazov-Larkin cor- 
rection on Fig. 2 (b). After transforming the Matsubara 
sum into a contour integral and going to dimensionless 
variables, it readsB: 
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where both the energy and the momentum cut-off have 
been set to unity. When renormalizing, we will followa: 
first integrate out a thin outer shell in the momentum 
space between 1 and 1 — 1 /b. Then rescale the momentum 
(q — » qb) to restore the cut-off, then rescale the energy 
(2 — > zb 2 ), the mass term and the temperature (T — > 
Tb 2 ) and, finally, integrate out an energy shell to restore 
the energy cut-off. At each step the quartic interaction, 
induces corrections to the mass term and to itself (seecl 
for details of the renormalization group equations). 

As a result, one arrives at the following transformation 
law for the Aslamazov-Larkin correction: 

Act[J] = b 2 - D Aa[J']+lnbf[J], 

where J denotes the mass term, the temperature and the 
quartic coefficient, J' denotes their renormalized values 
and D is the dimensionality of the sample. The precise 
form of f[J] can be easily obtained using the above renor- 
malization procedure. However, only two features of f[J] 
are important: (a) as long as the running value of T is 
smaller than the cut-off, f[J] has rather weak dependence 
on its arguments; this corresponds to the quantum renor- 
malization region; (b) in the classical renormalization re- 
gion, when T exceeds the cut-off, /[J] is proportional to 
T . Thus Act can be written as 



Act = 



In b* 



dxf[J{e*)]eV- D >> 



(4) 



where b* is the value of the rescaling factor at which 
the mass term reaches the cut-off and the scaling process 
stops. It is of the or der of T~ 3 / 4 in three dimensions 
and of the order of ^/ln (1/T)/T in twotl To evaluate 
(4), one also needs the value of b\ such that T(pi) ~ 1, 
which is bi ~ T^ 1 / 2 independently of the dimensionality. 
With these prerequisites, the answers in Table I follow 
as soon as one neglects the dependence of f[J(e x )} on all 
the couplings except temperature: 

f[J(e*)] = f[Te 2 % 

In three dimensions, the same result can be obtained 
just by replacing the "bare" T 2 mass term in (2) by its 
renormalized value T 3 / 2 , and then calculating the con- 
ductivity correction as per (3). 

Scaling analysis also allows to show that the vertex cor- 
rections are negligible. Their inclusion reduces to putting 
into the diagram of Fig. 2 (b) extra bubbles such as one 
on Fig. 2 (c). Each bubble contributes two Green's func- 
tions plus an integral over the energy and momentum. 
After an infinitesimal scaling transformation, this gives 
an extra factor of b 2 ~ D . Thus, in three dimensions, ver- 
tex corrections are irrelevant. In two dimensions, they 
appear marginal yet do not introduce any extra correc- 
tions. This can be established by using the Ward identity 
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and writing the renormalization equation directly for the 
current vertex ( see Fig. 3), and then inserting the solution 
into the corresponding expression for the conductivity. 

We would like to discuss theoretical and experimeH,- 
tal implications of the results. As envisaged by HertzH, 
under quite general circumstances the superconducting 
quantum critical point falls into the z = 2 universality 
class. Although in principle different behavior cannot be 
ruled out, it appears to require additional tuning of pa- 
rameters, as well as rather unusual features of the pairing 
phenomenon iteelf, such as the gap vanishing at the entire 
Fermi surface.Q 

In a weakly interacting disordered metal, To in (2) is 
of the order of the impurity scattering rate. In this case 
anomalous corrections due to quantum criticality in the 
Landau-Ginzburg regime are of the order of the weak lo- 
calization corrections S3 Therefore, in two dimensions our 
results are fully consistent only as long as all corrections 
are small and additive. Moreover, in two dimensions, the 
electron-electron interaction is known to lead to the lin- 
ear |tepperature dependence of the quasiparticle decay 
rateBH, regardless of closeness to the quantum critical 
point. Thus, for a weakly interacting two-dimensional 
system, the entire normal region corresponds to the two- 
dimensional disordered metallic regime (see Fig.l) with 
the quasiparticle decay rate proportional to the quasipar- 
ticle energy. 

However, we expect that in a strongly correlated sys- 
tem with suppressed weak localization effects, supercon- 
ducting quantum critical point still falls into the z = 2 
universality class and Table I may describe the actual 
state of affairs for all dimensionalities. In this case, Fermi 
liquid turns marginal only at T c = and crosses over to 
the normal Fermi liquid behavior as the system is driven 
away from the quantum critical point into the metallic 
phase. 

Finally, we should like to comment on the diagnostic 
opportunities, furnished by measurements at a supercon- 
ducting quantum critical point. In light of the discussion 
in the beginning of the paper, one is led to conclude that 
in a clean time reversal invariant system, observation of 
singular behavior at the superconducting quantum criti- 
cal point would mark a very peculiar phenomenon, as in a 
clean BCS superconductor suppression of T c completely 
eliminates pair fluctuations. Since any sample contains 
impurities, in reality the above conclusion refers to tem- 
peratures above the elastic scattering rate: T > t _1 . 
The latter can be extracted from the residual resistivity 
measurements of the sample. 

To summarize, we studied the properties of a su- 
perconductor near a quantum critical point driven by 
pair-breaking disorder. Superconducting fluctuations are 
characterized by a dynamical critical exponent z = 2. In 
a BCS superconductor at experimentally accessible low 
temperatures, the singular contribution to the conductiv- 
ity is proportional to y/T in three dimensions, but to 1/T 
in two dimensions. In a superconductor with the char- 
acteristic energy scale of the order of ef at the quantum 




FIG. 2. (a) Graph for the quasiparticle self-energy due 
to scattering by pair fluctuations. Solid lines denote the pair 
propagator; thin line denotes the quasiparticle Green func- 
tion, (b-c) Graphs for the conductivity correction: (b) the 
Aslamazov-Larkin graph (c) an example of a graph with ver- 
tex corrections. Scaling dimension of all such graphs is an 
integer multiple of (2 — D). 




FIG. 3. (a) Graph for the conductivity, including the 
current vertex, denoted by the dark triangle, (b) contribution 
to the free energy in the second order in the quartic term, 
which gives rise to renormalization of the current vertex. 

critical point, the contribution to the conductivity varies 
as T 1 ' 4 in three dimensions, and as In (1/T) in two. 
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